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(Dated: October 31, 2008) 

We present a study of the polarization of the T(15) and T(25') states using a 1.3 fb"'^ data sample 
collected by the DO experiment in 2002-2006 during Run II of the Fermilab Tevatron Collider. 
We measure the polarization parameter a = (ctt — 2aL)/{oT + 2aL), where ot and ol are the 
transversely and longitudinally polarized components of the production cross section, as a function 
of the transverse momentum (py) for the T(15') and T(2S). Significant pj-depcndcnt longitudinal 
polarization is observed for the T(IS'). A comparison with theoretical models is presented. 
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PACS numbers: 14.40.Nd, 13.88.+e, 13.20.Gd 



The production of heavy quarks and quarkonium states 
at high energies is under intense experimental and the- 
oretical study [H- The non-relativistic QCD (NRQCD) 
factorization approach has been developed to describe 
the inclusive production and decay of quarkonia 0] in- 
cluding high transverse momentum {px) S'-wave charmo- 
nium production at the Fermilab Tevatron Collider 
The theory introduces several nonperturbative color- 
octet matrix elements (MEs). These MEs are univer- 
sal and are fitted to data of the Fermilab Tevatron Col- 
lider 0] . The universality of the MEs has been tested in 
various experimental situations Q . A remarkable predic- 
tion of the NRQCD approach is that the S'-wave quarko- 
nium produced in the pp collision should be transversely 
polarized at sufficiently large pT @|. This prediction 
is based on the dominance of gluon fragmentation in 
quarkonium production at large px and on the ap- 
proximate heavy-quark spin symmetry of NRQCD 
Measurements of the polarization of prompt J/ip by the 
CDF Collaboration do not confirm this prediction [7[ . 

A convenient measure of the polarization is the variable 



a ((TT - 2(Tl)/{'7t + 2(Tl), 



(1) 



where ar and ul are the transversely and longitudinally 
polarized components of the production cross section. 
If we consider the decays of quarkonium to a charged 
lepton-antilepton pair, then the angular distribution is 
given by 



dN 



c?(cos ( 



(X 1 -I- a cos 



(2) 



where 6* is the angle of the positive lepton in the quarko- 
nium center-of-mass frame with respect to the momen- 
tum of the decaying particle in the laboratory frame. 

Quantitative calculations of the polarization for inclu- 
sive T{nS) mesons are carried out [§] by using the ME 
for direct bottomonium production determined from an 
analysis of Tevatron data 0]. They predict that the 
transverse polarization of T(IS') should dominate and in- 
crease steadily with p^ for p^ > 10 GeV/c and that the 
T(25') and T(3S') should be even more strongly trans- 
versely polarized. The fct-factorization model [lO| , using 
a semi-hard approach, predicts a longitudinal polariza- 
tion of T(IS') at p'^ > 5GeV/c [llj. In this context, 
the experimental measurement of the T polarization is a 
crucial test of two theoretical approaches to parton dy- 
namics in QCD. 

The DO detector is described in detail elsewhere 



The main elements relevent to this analysis are a central- 
tracking system, consisting of a silicon microstrip tracker 
(SMT) and a central fiber tracker (CFT), and muon de- 
tector systems. 



The data set used for this analysis includes approxi- 
mately 1.3 fb~^ of integrated luminosity collected by the 
DO detector between April 2002 and the end of 2006. We 
selected events where the T{nS) decayed into two muons. 
Muons were required to have hits in three muon layers, 
to have an associated track in the central tracking system 
with hits in both the SMT and CFT, and to have trans- 
verse momentum p!^ > 3.5GeV/c. In this analysis only 
events that passed a dimuon trigger, which requires two 
opposite charge muon candidates, were included in the 
final sample. We observed about 260,000 T{nS) with ra- 
pidity 1 2/^ |< 1.8 when fitting the dimuon invariant mass 
distribution as described below. 

Monte Carlo (MC) samples for unpolarized T(IS') and 
T(2S') inclusive production were generated using the 
PYTHIA [3] event generator and then passed through 
a GEANT-based [lj| simulation of the DO detector. The 
simulated events were then required to satisfy the same 
selection criteria as the data sample including a detailed 
simulation of all aspects of the trigger requirements. 

We fitted the dimuon invariant mass distribution in 
several intervals of p^ for a set of | cos 6**1 bins. A pre- 
vious measurement of the T(IS') cross-section by the DO 
experiment [l3| showed that a double Gaussian function 
is required to model the mass distribution of the T(IS') 
candidates. Studies performed on the T(IS') Monte Carlo 
sample suggest that a more sophisticated parameteriza- 
tion of the invariant mass distribution for some | cos 0* \ 
bins, where we observe non-Gaussian tails, is required. 
Two different parameterizations of the mass distribution 
were used, referred to as "data-driven" and "MC-driven" 
functions. The data-driven function has the advantage 
that no assumptions are made about how well the MC 
reproduces the true resolution. It consists of a double 
Gaussian function with equal means. The mean, widths, 
and relative fraction are free parameters. In contrast, the 
MC-driven function allows for a test of the effect of non- 
Gaussian components to the resolution that are observ- 
able in MC but are hidden in data by the detector reso- 
lution and the combinatoric background. Non-Gaussian 
tails arc implemented via a third Gaussian component 
with a floating mean to account for an asymmetric tail 
in the reconstructed T{nS) mass. The width and rela- 
tive fraction are taken from Monte Carlo. Figure [1] shows 
an example of a fit to the mass distribution for a single 
p^ and I cos 0*1 bin ignoring or including non-Gaussian 
tails. The signal consists of three mass peaks, the T(IS'), 
T(2S'), and T(3S') where the mass differences were fixed 
to the measured values [11]. The background was mod- 
eled with a convolution of an exponential and a polyno- 
mial function. The degree of the polynomial was chosen 
to be between one and six depending on the complex- 
ity of the shape of the background. The values in 
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FIG. 3: Comparison of data (points) and Monte Carlo (solid 
histogram) for T(IS') with 2 < < 4GeV/c: a) momentum 
of T(IS'), b) polar angle of T(15'), c) angle between muons, 
d) |cos6l*|. 



FIG. 1: [Color online] Signal extraction from the dimuon in- 
variant mass distribution for events in the 0.4 < | cos 9* \ < 0.5 
region. a,c) 2 < < 4GeV/c; b,d) 10 < < 15GeV/c. 
Dashed curves are the combinatoric background. 



'0.2 



^0.1 



~ "■ - 1 


a) 







0.2 



0.1 



b) 



0.2 0.4 0.6 0.8 .1 
IcosQ I 







0.2 0.4 0.6 0.8 .1 
loose I 



FIG. 2: [Color online] Monte Carlo | cos9*\ distributions after 
all selection requirements for different a values: —1 (dashed 
histogram), (solid histogram) and +1 (dotted histogram), 
a) < p? < 1 GeV/c, b) p? >15GeV/c. 



Fig. [T] do not allow us to differentiate between the two 
approximations and hence we average them. 

The data were divided into bins in and | cos 6* \ . For 
each of these bins the numbers of T(IS') and T(2S') can- 
didates were extracted from the mass distribution. The 
number of T{3S) candidates was insufficient to extract 
angular distributions. 

Polarization was not taken into account in the Monte 
Carlo generation. To compare them with data we calcu- 
lated for each event the weight Wa, which converts the 
initial Monte Carlo | cos 9* \ distribution with a = to 
a distribution with the chosen a. Figure O shows the 
sensivity of the DO detector to the T(IS') polarization 
for the lowest and highest p^^"^'' intervals. The pythia 
simulation does not accurately model the kinematic dis- 
tributions of T{nS) production at the Tevatron (e.g., the 
p^'-"'^-' distribution). To correct the Monte Carlo dis- 
tributions, we introduced additional weights to improve 
the agreement with data of the T{nS) momentum dis- 



tribution. Instead of the weight Wa in our algorithm, we 
used the weight w = WaWprWpr, where WpT and Wpr are 
weights to achieve agreement between data and Monte 
Carlo distributions of and . After this reweighting 
procedure, we obtained good agreement between data 
and MC for the T(nS) and muon kinematic distribu- 
tions. An example for Y(1S) with 2 < p^ < 4GeV/c, 
using the MC-driven fit, is presented in Fig[3l All data 
distributions were derived by estimating the number of 
T(IS') events from a fit to the dimuon mass distribution 
for the corresponding bin of the histogram. 

The systematic uncertainties on a for T(IS') are sum- 
marized in Table HI Values of a were found for several 
p^ intervals, using both parameterizations (data-driven 
and MC-driven) of the dimuon invariant mass distribu- 
tion for the signal. Both a measurements are averaged 
and one half of the difference between them is assigned 
as systematic uncertainty due to the signal model. The 
uncertainty in the background was estimated by varying 
the mass range of the fit and the degree of the polynomial 
used to parameterize the background. The MC simula- 
tion does not reproduce exactly the mass of the T(IS') 
peak, which differs by about 40MeV/c^ from the PDG 
value. The effect on the a determination was estimated 
and shown in Table [J under "muon momentum." Fi- 
nally, the systematic uncertainty due to the trigger sim- 
ulation has also been considered and shown in Table [H 
The T(IS') polarization was calculated assuming that it 
is constant within a given p^ bin. This assumption leads 
to a small bias in the measured a that is estimated by 
reweighting the simulation using the observed p^ depen- 
dence of a. The final measured a is corrected by a factor 
ranging between —0.03 and -1-0.06, depending on pj. 

Figure |4] shows the measured a as a function of pj, 
for T(IS'). Note that the bin for 14-20 GeV is not sta- 
tistically independent from the adjacent bins. The ar- 
row indicates that the highest p^ interval considered. 
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TABLE I: Systematic uncertainties on q for T(IS'). 



Source 


Uncertainty on q" 


pI ' [GeV/c] 


Signal model 


0.01-0.15 


1-2 


Background model 


0.04 - 0.21 


- 1 


Muon momentum 


0.00 - 0.06 


0-1 


Trigger simulation 


0.00 - 0.06 


>15 



"For all intervals 

''Interval with maximal uncertainty 



D0, 1 .3 fb"'' 




12.5 15 17.5 20 
Pj of r(1S) [GeV/c] 

FIG. 4: [Color online] Dependence of a on for inclusive 
T(1S) candidates. Black circles are data. The band is the 
NRQCD prediction @]. Curves are two limiting cases (see 
text) of the fct-factorization model pT|. 



> 15GeV/c, does not have an upper limit. The un- 
certainties are the systematic and statistical uncertain- 
ties added in quadrature. Also shown are the NRQCD 
prediction Q (yellow band), and the two limits of the 
fc(-factorization model [ll| (curves). The lower line cor- 
responds to the quark-spin conservation hypothesis, and 
the upper one to the full quark-spin depolarization hy- 
pothesis. The previous measurement by CDF of the po- 
larization of T(IS') with rapidity | |< 0.4 is consistent 
with a equal to zero We expect the CDF and DO 

results to be similar and we have no explanation for the 
observed difference. We also extracted the polarization 
of the T(2S), which is shown in Fig. [5] along with the 
NRQCD predictions Values of a for statistically in- 
dependent intervals, shown in Fig. [5] and Fig. O are 
given in Table HIl 

In conclusion, we have presented measurements of the 



TABLE II: Measurements of a for T(15') and T(15'). 



Pt [GeV/c] 



«[T(lg)] 



"[T(2g)] 



0- 1 

1- 2 

2- 4 
4-7 
7-10 
10 - 15 

>15 



0.04 ±0.27 
-0.41 ± 0.20 
-0.54 ±0.17 
-0.55 ±0.10 
-0.45 ±0.14 
-0.34 ±0.14 

0.25 ±0.19 



-0.04 ±0.54 
0.28 ± 0.37 
0.04 ±0.18 

-0.37 ±0.21 
0.09 ±0.32 
0.32 ±0.27 
0.55 ±0.58 



B 1 F T 

° D0,^ .3 fb"^ 

0.8 : 

0.6 : 

0.4 : 

0.2 

i| 
-0.2 
-0.4 : 
-0.6 : 

"°'^0 2.5 5 7.5 10 12.5 15 17.5 20 

Pt- of r(2S) [GeV/c] 

FIG. 5: [Color online] Dependence of a on for inclusive 
T(25') production. Circles are our data. The band is the 
NRQCD prediction [|. 



polarization of the T(15) and T{2S) as functions of 
from OGeV/c to 20GeV/c. Significant pT-dependent 
longitudinal polarization is observed for the T(IS') in- 
consistent with NRQCD predictions. At >7GeV/c 
the fraction of transversely polarized T{2S) particles is 
higher than in T(IS') at the same value of p^, in agree- 
ment with NRQCD predictions. 
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